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We have constructed the bag model having a central constant color field. The motion 
of the quark is studied in this bag and the Dirac equation is solved for it. The energy 
O spectrum found has a branching due to the interaction of the quarks with the color 

background. It is pointed out that this model can be applied for taking into account, 
in mass spectrum of the hadrons, the coupling of the constituent quarks with the gluon 



condensation as the interaction with the color background. 



Introduction 



There are a constituent quark model and some bag models describing hadron states by means 
of the dynamics and kinematics of the constituent quarks [1-4]. These models take into 
account the gluon and quark condensation existing in the QCD vacuum [5] as a medium in 
which the constituent quarks move and the interaction of these quarks with condensates are 
reduced acquiring additional mass by them. The vacuum average values of the chromofield 
components (S a £ a ) and (TL a Ti a ) of the gluon condensate in the QCD vacuum were estimated 
in [5-7]. These estimations are used for the study of vacuum energy and structure of the 
bag models, which connected with the bag constant [8,9]. Here we aim to consider the 
gluon condensate as a color background field and take into account the coupling of the 
constituent quarks with it as an interaction with the external color filed. Naturally, such a 
coupling will rise in the constituent quarks energy spectrum which is not consisting in only 
acquiring additional mass term and will result in a splitting of this spectrum. In order to 
provide the rotational invarince we propose that spatial components of this background are 
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equal and constant for simplicity. The background field thus defined will play the role of 
a central constant color force. In order to realize our proposals we shall construct a model 
like Bogolubov's bag model [3], but having a constant color field in the center. According to 
this model the colored constituent spinor particle moves in the constant central color field 
and its motion is limited by an infinite spherical well. The condition of an infinite spherical 
well is imposed in order to ensure the confinement property of the colored constituents. 
The Dirac equation should be solved for this model under this boundary condition on wave 
functions of the constituent quark. For this aim we can use the constant non-abelian vector 
potentials found in [10] and the approach used for solving the Dirac equation in the field 
given by such a kind of potentials [11]. This will give us the distribution of constituent 
quarks inside the bag and their energy spectrum, which will contain the contribution of its 
color and spin interactions with the background field. Finally, the components of external 
field can be identified with the ones of the gluon condensate and the estimations for the 
color components of the gluon condensation (S a S a ) and (7i a 7i a ) can be setted in the energy 
spectrum found for the quark. This will make this spectrum useful for a comparison with 
the spectra of hadron states. 



1 The Dirac equation 

For a study of the problems in a constant non-abelian background field it is convenient to 
introduce the constant non-commuting vector potentials 1 . Let us choose the vector potentials 
in the framework of the SU C (3) color symmetry group in the following way: 

4 a) = ^n, Af ] = ^5 ja (j = 1, 2, 3) for a = 1, 2, 3 

4«)=0 for a = 4,5,6,7,8, 1 ; 

where r, T\ are constants and 5j a is the Kronecker symbol. The corresponding field strength 
tensor Fffl = gf abc A^A^ has the following chromoelectric (£f^) and chromomagnetic 

(Hj 1 ^) components: 

£V) = 0, = g^TTi, 4 3 ) = -g^FTu W« = gr, Uf = 0, = 0, 
£$> = -Qy/rn, = 0, £® = 9y/rn\ = 0, Uf = gr, Hf = 0, (2) 
= 9^ £J 2) = -9 y/fn, Sf ] = 0; HP = 0, Hf = 0, Hf = gr. 

Here g is the color interaction constant. All other color components of Fj$ are zero for 
a = 4, 5, 6, 7, 8. From (2) we see that the field (1) has equal magnitude of spatial components 
£2 = Efsf = 2g 2 rr 1 and U) = Jifuf = g 2 r 2 , which are constant as well. So, in 
ordinary space the strength vectors of the chromomagnetic and chromoelectric fields are 2 : 

7i = V3grlt, = gV&rri it, (3) 

where rf is the unit radius vector in ordinary space. 
1 See [11] and references therein. 

2 We assume the chromoelectric field has negative projections in color space. 
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The Dirac equation for a colored particle minimally coupled with the external color field 
(1) can be written as follows: 

(7^ - M) * = 0, (4) 

where F M = p^ + gA^ = p^ + gAJ^^y , and the A a are Gell-Mann matrices describing the 
particle's color spin. Equation (4) can be written in terms of the Majorana spinors (ft and x 



X 



in the more suitable form for us: 



(Vp,) 2 v>= (p 2 -m 2 )^ 



(5) 



where the Pauli matrices Oi describe a particle's spin. Here and afterwards tp means (ft or 



X- The two spin components of the Majorana spinors ip 



1p2 



transform under the 



fundamental representation of color group SU C (3). That means that each spin component of 
the wave function ^ ij2 has three color components describing the color states of a particle: 



^1,2 = 



Writing down the expressions of F M and Affl in equation (5) we get an explicit form: 



> 2 ' - 2 ' 3 ^ +gT h a p a -^a a X a )^=(E 2 -grh a E + ^ 



~p> + M + 



V, (6) 



where i? is the energy of particle 3 . Equation (6) turns into following system of differential 
equations for the color components ift[ a 2 - 



A + gr-2V 3 - \g 2 r) ^ + gr\ (Vi - %V 2 ) 4? 
(A - gAv 3 + \g\) 4? + gr\ {V x + iV 2 ) ^ 

(y 2 + M 2 ) ^ 3) 



A 



gT*V 3 + \g 2 r) V4 



(1) 



g T -2 (V! - iv 2 ) n 

A - , r r- P, - y 2 r ) tjj {2) + gr\ (Vi + lP 2 ) 4 

(f 2 + M 2 ) i). 



(2) 



.(1) 



(3) 





g 2 Tip 2 

E 2 i,f ] 
g 2 Tip[ 


e 2 4 3) 



(i) 



(2) 



(7) 



where the operators A and Vj denote A = ~f 2 + M 2 + \g 2 (r - nj-E 2 , Vj = pj + ^E. The 
equations in the system (7) mix the different states ift[ a 2 and we need separeted equations for 



i Since the field (1) does not depend on time, the states are stationary. 
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each of these states. From the system (7) we find that equations for all states i/j[^ 2 (i — 1,2) 
have the same form: 



<ptl = 0, 



(8) 



which possesses rotational invariance. Since the operators A and Vi commute, the operator 
in first square bracket commutes with the second one. This allows us to divide the (8) into 
two equations, both keeping rotational invariance: 



2 \ 2 

9 t\ 



9 2 rV 



Vfi = o, 



(9) 



A + 



-g 2 r[V +g 2 r 



1$ = o. 



(10) 



Equations (9) and (10) can be solved separetaly and the set of solutions (8) will consist of 
the solutions of (9) and (10). Let us consider equation (9). Acting on this equation by the 

operator ~fi we get the same equation for the function £^ — ~^ ^1*2 as f° r V'i'L i- e - : 



This means the functions an d 1P12 differ only by a constant multiplier k' 2 : = k' 2 ipi :2 



'.(0 



& = 0. 



12 . t(0 



or 



= 



(ii) 



In other words, since the operator V commutes with the square bracket operator in (9), 
they have the same set of eigenfunctions The same claim is order for the operator ~p" 2 : 



2./.W 



(12) 



Thus, we can solve (12) instead of (9). Obviously, (12) keeps the rotational invariance 
property of equivalent equation (9) and so is easily solved in a spherical coordinate system 
using the separation ansatz [12]: 



(13) 



Here r = y/x 2 + y 2 + z 2 , I and m are the orbital angular momentum and chromomagnetic 
quantum numbers, 6, ip are the polar and azimuthal angles, respectively. The spherical 
functions Y t m (9, ip) are expressed by means of the Legendre polynomials Pj m ^ (cos 9) : 



Y i m (8, V>) 



{21+1) (l- 


m )! 


4tt (/+ 


m )! 



(cos 9) e 



imip 



(14) 
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and define the "s" , "p" , "cf ' , "/" , • • • orbitals well-known in quantum mechanics. The 
equation for the radial part R(r) is as used in many quantum mechanical problems possessing 
rotational invariance [12,11]: 

|i*M + + (* 2 " ^) «M = o- < 15 > 

With the notations Q (r) = y/rR(r) (15) turns into Bessel's equation for Q ( r ) : 

Q" (r) + V (r) + ^ 2 - j Q (r) = 0. (16) 

The function i?(r) must be finite on r — > 0. This means that for the solution of (15) we 
should choose the Bessel function of the first kind: 

R i (r) = -^=J l+1/2 (kr). (17) 

Thus, we conclude that states obeying (9) are the "s" , "p" , "cT , "/" , • • • orbitals corre- 
sponding to the different values of the quantum numbers / and m and the motion of the 
constituent quarks in any color and spin state takes place on these orbitals, which are the 
same for all these states. In other words, the angle distribution of quarks inside bag is the 
same for any color and spin state and the same, for instance, as ones of electrons in atoms 4 . 
The radial distribution of these quarks is the same as ones of freely moving particles enclosed 
in a sphere [12]. Since, the external field (1) does not depend on r, we have obtained the 
same expression for the solutions \ (~^) as f° r a freely moving particle enclosed in a sphere, 
differing only by the expression of the k 2 constant. 



2 The energy spectrum 



Using (3) the operator V can be written in the following form: 

E 



V 



Then the action of 7^ operator will be 



2,/,(0 



IE 



VQg' 



t 



cosa H E w\ o, 

T 



where 



^■01*2 means the directional derivative [13]: 



\\dxj \dy J \dz / 



(18) 



4 Since constituents in this model are considered to be non- interacting with each other, for the two-particle 
case we have the angle distribution of quarks in mesons. This distribution will have the same shape as ones 
in two-electron atom. In this sense mesons like hydrogen atom, differing only by the radial distribution in 
the constant field approximation. 
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Here a is the angle between the momentum vector V "0 and the chromoelectric field vector 
'. Since the particle periodically moves on orbitals the angle a varies in symmetric limits. 
That means the average value of cos a during one period is zero. Actually, the term, which 
contain cos a in (18) is proportional to the work done by chromoelectric field on the particle 
in its motion in this field. It is easily seen the net work of this field during one period is zero, 
while its momentary value is not zero 5 . So, if we average (9) and (10) over the time during 
one period, the term proportional to cos a will drop out. We shall find the average value of 
energy spectrum 6 . Therefore, the average value of (18) during one period is equal to: 

= -V 2 ^ 2 + T ±E^% = k*4% + = k'*4% (19) 

If we take (19) and (12) into account, we obtain from (9) the following equation for the 
constant k 2 : 

(k 2 f + 2k 2 [m 2 - E 2 + \g 2 (r - 3n) - \g 2 r^ + [m 2 - E 2 + \g 2 (r - Zntf -g 2 Tl E 2 = 0, 



from which one finds the relation between the constant k 2 and the energy spectrum E~ 

2 



2 . 




Now we should impose the boundary condition meaning the confinement property of 
constituent quark on its wave function. Since (12) is the Laplace equation we can impose 
the Dirichlet boundary condition if)[^ 2 (r = r ) = or Ri(r ) = 0, which means that we 
enclose the particle motion by a sphere with a radius r . Here r agrees with the half of 
hadron size. This boundary condition establishes the following relation between the values 
of k and the zeros of the Bessel function J;+i/2 (kr) : 

kr = al N) (21) 

which means the quantization of the k values. According to the relation (20) the energy 
spectrum of particle is quantized as well. Plugging (21) in (20), we find the first two branches 
of the quantized energy levels of the spectrum: 



2 

2 



1 




+ M 2 + - A g 2 (t - 2n) =F ^ 



T / OLi 



J 



(22) 



In (22) and (21) iV labels the sequence of zeros of the Bessel function iV = 1,2,3 and is 

called the radial quantum number [12]. So, the value of angular momentum I determines the 
series of the energy spectrum of the particle and the radial quantum number iV determines 
the energy levels in this series. Thus, the finiteness condition imposed on the motion of 

5 We observe the same situation as in an electron's motion in the field of nucleus. 

6 For "s" orbitals cosa=0, since the momentum and chromofield vectors are perpendicular on every 
moment in time. So, the momentary value of energy coincide with its average value for these orbitals. 
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the particle, because of its confinement property, of course, leads to the quantization of the 
energy spectrum. 

As is seen from (21) the constant k gets the same values for the various branches (k 2 ) 1 2 : 



k. 



(AO 
(AO = Oj 

?~0 



The radius a of the turning point of the particle could be found from (15) using the maximum 
condition on the radial function R'(a) = and is equal to: 



(N) 

a) = r 



1(1 + 1) 

a) 



These radia, as the energy spectrum and k, are quantized and are determined by the quantum 
numbers / and N and do not depend on the field intensities. So, there is no difference in 
values of the radia with the case of motion in a pure chromo magnetic field [11]. 

We have solved (9) and found the corresponding two branches of the energy spectrum. 
In the same manner we can solve (10) too. For this equation the function denotes 

d*2 — (j^ + 9 2r ^j "01*2 an d the following equations are equivalent to (10): 

(V 2 + 9 2 r) Vfi = K^%, (V 2 + g 2 r) ^ = (23) 

The relation between the expressions of k and K is obvious K 2 = k 2 + g 2 r, K' 2 = k' 2 + g 2 r 
and the relation between K and K' is analogous to (19): 

K' 2 = K 2 + -E 2 . (24) 

T 

Taking (23) and (24) into account, we get from (10) an algebraic relation between k 2 and 
E 2 : 



(k 2 + M 2 + \g 2 (r - Tl ) - E 2 + - g 2 r (k 2 + ^E 2 + g 2 r) = 0. 



(25) 



Of course, the solution of (23) is (13) with (14) and (17). The Dirichlet boundary condition 
Ri( r o) — applies to (10) too. In the result of the quantization (21), from (25) we find the 
other two branches of the energy spectrum: 



(EDI = J (1 + i) ((4?) 2 + 9 2 -) + + b 2 (r - 2n) =F \grj 



ti \ \ ro 



+ g 2 r 



(26) 



As seen from (22) and (26), the energy spectrum of the quark contains the contribution of 
quark's interaction with the background field and this energy spectrum is quantized due to 
the finiteness of motion of quark. Another difference with the existing bag models [1-4] is 
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the branching of the energy spectrum, which is the result of the color and spin interaction of 
the particle with the background field and is not determined by these quantum numbers. We 
have four color and spin states \ °f the quark and four energy branches (£7^)1,2,3,4, but 
we do not have one-to-one correspondence of these states and the branches. All these states 
could get energy from the any of these branches. This is the difference between the branching 
of spectrum and the splitting levels. In fact, here we have obtained the branching of energy 
spectrum instead of splitting levels. Such a branching takes place in a pure chromomagnetic 
field case [11,14,15] and is connected with the existence of conserved operator, which contains 
the color spin 4j- and the spin ^ operators. There is no need to introduce the color states 
in bag models, which do not deal with the color interaction of particle with the background, 
and so, this branching of energy spectrum does not occur in those models. Note that the 
energy levels in these branches are determined only by the quantum numbers / and N, and 
so they are 2/ + 1 fold degenerate in the quantum number m, which occurs on motion in any 
central field. 

The states ip^ correspond to the states of a colorless particle, and (12) for them (k 2 = E 2 — M 2 ) 
has the solution (13) with (15) and (17) too. The energy spectrum of these states is: 




As an application of this model we can propose the central color field is the background 
field of gluon condensation in QCD vacuum. Then we can use the estimation for vacuum 
average values of chromoelectric and chromomagnetic field intensities of this condensation 
[6,7]: 

(0 I g 2 t a t a I 0) ~ - (700 Mevf , (0 | g 2 lTlT | 0) ~ (700 Mevf . 

Identifying these estimates with the field intensities in our model g 2 E = Qg A TTi and g 2 H = 
Sg^r 2 , we see that constant T\ should be taken equal to T\ = —\t. This enable us to evaluate 
the g 2 r constant g 2 r = ^ (700 Mev) 2 . Then the energy spectra (22) and (26) simplify as 
follows: 




In order to estimate quark's energy using the spectra (28) and (27) we can set in them 
ro = R — 1, 7 fm [3] or another estimation for hadron's size. Since constituent quarks in the 
framework of this model are considered non-interacting with each other, the energy of any 
two- or three-quark system will be the sum of the energy of separate quarks. As an exam- 
ination of this model, the total energy of the bag can be compared with the corresponding 
hadron state. It should be noted that the selection rule Am = 0, ±1 and Al = ±1 for quan- 
tum transitions between the energy levels, which exists in central field problems [12], occurs 
for this problem as well and is useful for a calculation of the energy emitted by the excited 
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bag. There is not preferred rule on the energy branches for this calculation. A comparison 
of the energies emitted by bag and by the hadron could serve as a another verification of the 
role of the gluon condensate as constant central color field. 
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